Non-uniqueness of phase shift in central scattering due to monodromy 
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Scattering at a central potential is completely characterized by the phase shifts which are the 
differences in phase between outgoing scattered and unscattered partial waves. In this letter it is 
shown that, for 2D scattering at a repulsive central potential, the phase shift cannot be uniquely 
defined due to a topological obstruction which is similar to monodromy in bound systems. 
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Introduction. — In physics it is often crucial to find 
"suitable" coordinates. One important set of coordinates 
in Hamiltonian mechanics are the action-angle variables 
[l|. Action-angle variables played a crucial role in the 
development of early quantum mechanics in the Bohr- 
Sommerfeld quantization rule. The existence of action- 
angle variables was already addressed by Einstein in 1917 
[2| . He argued that a quantization of actions only works 
if the mechanical system is integrable, i.e. it has as many 
independent constants of motion in involution as degrees 
of freedom. The existence of action-variables was made 
more precise in the Liouville- Arnold theorem [3]: if a 
connected component of the common level set of the con- 
stants of motion is regular (i.e. the gradients of the con- 
stants of motions are everywhere linearly independent on 
the component) and bounded then it has the topology of 
a torus and action-angle variables exist in the neighbor- 
hood of the torus. The angles are the coordinates on 
the individual tori (which are just Cartesian products 
of circles) and the actions change from torus to torus 
in a smooth way. This theorem is local near a regular 
torus. The question of if and how these "local action- 
angle variables" fit together globally has been ignored for 
a long time. Recently it has been shown that there can 
exist topological obstructions to the global uniqueness of 
action-angle variables [4]. Consider a family of regular 
invariant tori in phase space that starts and ends with 
the same torus. This family may be a non-trivial torus 
bundle. Similarly, a Mobius strip is a non-trivial inter- 
val bundle: even though every local piece of it is just a 
rectangle, globally it is twisted. If the action variables 
are changed by a non-trivial uni-modular transformation 
after one (mono) circuit (dromos) through a family of 
regular invariant tori, then the system has monodromy. 
Monodromy implies that the action-angle variables do 
not give global coordinates. It is a topological obstruc- 
tion, because the twist in the bundle cannot be removed 
by smooth deformations of the bundle. By contrast, if 
the loop of regular tori can be contracted (passing only 
through regular tori) it cannot have monodromy. In this 



way monodromy is related to a non-regular level set of 
the constants of motion that is not a torus. Such critical 
sets appear in phase space where the gradients of the con- 
stants of motion are linearly dependent. The most promi- 
nent example of a critical set that causes monodromy is 
a pinched torus 5]. It exists in integrable systems with 
two degrees of freedom that have an unstable equilibrium 
point of focus-focus type, i.e. with eigenvalues of the form 
a ± iu>, —a ± ico, a ^ 0. One of the simplest examples 
with this type of monodromy is the spherical pendulum 
[5( . The quantum version of this phenomenon [6|] explains 
why there is no global quantum number assignment for 
the hydrogen atom in external fields [7|, the Hj molec- 
ular ion [8(, the rovibrational spectrum of CO2 [9], and 
other systems [10j . 

In this letter we study the implication of the un- 
bounded analog of monodromy in scattering problems. 
Elastic scattering at a central potential is completely 
characterized by the phase shifts which are the differ- 
ences in phase of outgoing scattered partial waves and 
outgoing unscattered partial waves. For a planar system, 
a partial wave (x,y\l,p) with angular momentum I and 
asymptotic momentum p at infinity gains a phase 5(1, p). 
In fact, the action of the scattering matrix Sona partial 
wave is 



S\l,p)=exp(2i5(l,p))\l, P ) 



(1) 



All physical quantities such as scattering cross sections 
and amplitudes can be expressed in terms of 5(1, p). 

The phase shift is positive for attractive potentials and 
negative for repulsive potentials. If there is no interac- 
tion then 5(1, p) = 0. Similarly, the phase shift vanishes 
in the limiting case of large p where the potential can 
be ignored due to the dominating kinetic energy. The 
common procedure to define the phase shift therefore is 
to smoothly continue 5(1, p) from large to small p. We 
will show that for smooth repulsive potentials, there is a 
topological obstruction to this procedure, and as a con- 
sequence, the phase shift cannot be uniquely defined. 




FIG. 1: Phase portrait (r,p r ) withp 2 = 2fiE-l 2 /r 2 -2fiV(r) with 
E = 3, I = 1, and V(r) = (outer curve) and V(r) = a/(l + (br) 2 ) 
with a = 20 and 6 = 1 (inner curve). The shaded area is equal to 
AW denned in J2}. 



Nonuniqueness of phase shift. — Consider a smooth 
repulsive central potential V(r) with V(r) — > suffi- 
ciently fast for r — ► oo. The Taylor expansion at the ori- 
gin is V(r) = E c -(ia 2 r 2 /2+0(r' i ) with E c > 0. A semi- 
classical expression for 8(1, p) can be obtained from the 



WKB method [ll|. Assuming 8(1, p) — > for p — > oo the 



WKB approximation yields ^wkbCp) = AW(l,p)/(2H) 
where AW(l,p) is the difference of the radial actions with 
and without potential, 



AW(l,p) = W(l,p) - W'(l,p) :-- 



2 I yV - l 2 /r 2 ~ 2[iV(r) dr - 2 



(2) 
\J p 1 — I 2 /r 2 dr . 



Here r$ and r' are the classical turning points with and 
without potential, i.e., ro is the largest nonnegative root 
of r 2 p 2 — I 2 — 2[ir 2 V(r) or zero if I = in combination 
with p > p c = (2^iE c ) 1 / 2 , and r' — \l\/p. The difference 
AW is finite while the individual integrals diverge (see 
Fig. [J). 

Surprisingly the function AW(l,p) is not globally 
smooth: it is not differentiable at I — when p < p c . 
To illustrate this we show contours of AW in Fig. [2^. 
Consider the derivative of AW with respect to I, 



-AW(l,p) = — — 



-1 



z^Jzp 2 -I 2 - 2[izXJ(z) 



: dz — I 



Z\j zp 2 — I 2 



(3) 
dz. 



Here we substituted z — r 2 and let U(z) = V(y/z). In 
contrast to the integrals in ((2|) their derivatives with re- 
spect to I exist individually. The second integral is ele- 
mentary and gives — sgn(/)7r. 

For general I and p, the integral dW(l,p)/dl in ([3]) 
depends on the potential U. Interestingly, for / = 0, this 
is no longer the case. The limiting case I — > is tricky: 
depending on whether p > p c or p < p c the branch point 
zq of the square root in the integrand either collides or 
does not collide with the integrand's pole at z = as I — > 




FIG. 2: (£,p)-plane with contours of (a) AW defined in J2} and 
(b) AW defined in ((10). The bold dot marks (l,p) = (0,p c ). The 
potential is the same as in Fig. [T] 



0. The collision of the branch point and the pole leads 
to the divergence of the integral and the question arises 
of how this divergence is compensated by the vanishing 
of the prefactor I. At z — the argument of the square 
root has the Taylor expansion 



;2 i / 2 2\ , 2 2 2 

l + (p — p c )z + [i a z 



0(z 3 



(4) 



We note that in case a = higher order terms can 
be included and do not lead to a substantial change of 
the following argument. For p < p c , the linear term of 
this expression has a negative coefficient that does not 
depend on I. Hence, the collision of zq and zero as / — > 
does not take place, and the integral is not critical in this 
limiting case. Due to the prefactor I in (J3j> we thus have 
dW(l,p)/dl — ► as / — > and accordingly, the left and 
right hand derivatives of AW(l,p) with respect to / at 
zero are 



l-^OT ol 



T7T (p < p c ) ■ 



(5) 



When p > p c the coefficient of the linear term in (J4j) is 
positive and the collision of the branch point and the pole 
does take place. Consider the integral in the complex 
plane. For I ^ we define the integration path C as 
shown in Fig. [3] This gives 



dW(l,p) 
dl 



= 1 



-1 



c 2zy / zp 2 - I 2 - 2[izU(z) 



(6) 
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FIG. 3: Complex z-plane with integration paths for the differential 
—dz/(2z(zp 2 — I 2 — 2fizll(z)) 1 ' 2 ) which has a pole at the origin 
2 = (marked by the cross). The square root is real along the 
branch cut which extends from the turning point zq (marked by the 
dot) along the postive real axis; it is positive 'above' and negative 
'below' this branch cut. The integration path C (solid line) is 
equivalent to the composition C = C\ + C2 (dashed line) . 
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To study the limit I — > we deform the integration path 
C by wrapping it over the pole at zero and compensate 
the capture of the pole by adding a small closed integra- 
tion path that encircles the pole in opposite direction. 
Thus we decompose C — C\ + Ci with C\ and C2 as 
shown in Fig. [3J The integal ^ thus becomes 



dW(l,p) 
dl 



2 



= > I 



fe=i 



-1 



c k 2zsjzp 2 - I 2 - 2fizU{z) 



(7) 



The integration path C\ is not critical for I — > 
and due to the prefactor I the contribution to dW/dl 
vanishes for I — > 0. For the choice of the branch 
of the square root explained in Fig. [3] the differential 
-dz/(2z v / zp 2 - I 2 - 2^zU(z)) has residue -i/(2\l\) at 
2 = 0. The integral along Ci thus leads to the contribu- 
tion 



-1 



c 2 2zy/zp 2 -l 2 -2fizU(z) 



dz = sgn(/)7r . (8) 



We thus find that for p > p c , the contribution of 
dW{l,p)/dl and dW'(l,p)/dl to dAW(l,p)/dl cancel 
each other for I — ► and accordingly, 



Im^AW(i,p) = (p>p c 



(9) 



One might think of removing the kink in Fig. [5^ by 
'smoothing' AW(l,p) according to 



AW(l lP ) = 



AW(l,p) for/<0 
AVF(Z,p)-27r? for Z >0 



(10) 







This however introduces a kink at the segment of / 
where p > p c , see Fig. [2Jd. 

The WKB method gives 2h8 = AW where / = mH, 
m G Z. The values of 5 are only relevant mod7r, see 
(fT]). In other words, the function exp(2i<5) is (locally) 
periodic. In order to study this periodicity we consider 
the values p — kh such that S(mh, kh) = mod 7r, see 
Fig. [U The function exp(2i6) is not globally periodic 
because of the singularity at (to, k) = (0,p c /h). This can 
be seen by transporting a unit cell in the lattice around 



FIG. 4: Lattice of zeros (empty circles) of the phase shift <5 mod it 
in the plane m = l/h, k = p/h, and parallel transport of a lattice 
cell about the singularity (m, k) = (0,p c /K) (filled circle). The 
potential is the same as in Fig. [T] h = 0.25. 



the singularity. The lattice cell crosses the line to = 
according to the modified AW of (fTT)|) which is smooth 
for p < p c , while the original AW is smooth for p > p c . 
Thus in the presence of a repulsive localised potential 
the phase shift S cannot be globally defined. We call this 
phenomenon quantum scattering monodromy. 

Notice that the derivative of the phase shift with re- 
spect to the energy give the eigenvalues of the Wigner- 
Smith time delay matrix Q = —ihS~ 1 dS/dE. Semiclas- 
sically, this derivative is given by dAW(l,p)/dE = AT 
which is the classical time delay. This derivative is 
smooth everywhere apart from the point (I, E) = (0, E c ). 

The WKB approximation @ does not account for the 
collision of the classical turning point and the singularity 
of the effective potential when I — v 0. An asymptotic ex- 
pansion of the exact solution of the radial wave equation 
shows that nevertheless the error is less than 1%. 

Classical explanation. — The classical interpretation of 
dAW(l,p)/dl is the angle of deflection. Consider the 
polar angle if between the incoming and outgoing orbit 
in configuration space, 



<p-. 



f°° I 
tpdt= —zdt. 

J-00 Hr z 



Substituting di = dr/f gives 

<P = / -^-dr- 
Joe V T r 



00 I 1 , 
— r^-dr 

1 fir z r 



ro ry/2[iEr 2 -I 2 - 2fir 2 V(r) 



(11) 

(12) 
dr (13) 



where ro is the turning point. We used r < in the 
first integral and r > in the second integral in (JT2J) . 
Up to the sign the derivative dAW(l,p)/dl in ([3]) thus 
coincides with the angle of deflection Aip between the 
scattered and the corresponding unscattered orbit. 

Let us now follow the deflection angle Aip for pairs 
of angular momenta (or equivalently impact parameters) 




p + 



(c) 

p 



(b) 
(a) 




FIG. 5: Orbits in configuration space coming in from y = — oo for 
four different pairs of angular momenta and asymptotic momenta 
marked as bold points on the path in the (Z,p)-plane that encircles 
the critical point (l,p) = (0,p c ). 



and asymptotic momenta (i.e. energies) along a closed 
path in the (7,p)-plane that encircles the critical point 
(l,p) — (0,p c ), see Fig. [5] Consider orbits that come in 
from y — — oo. For I = and p < p c , the particle comes 
in along the y-axis and slows down. Due to insufficient 
energy it cannot overcome the potential barrier and so 
turns back towards y = — oo. This orbit has Aip = n. 
If / is increased to a value / = Z( a ) > (keeping p < p c 
fixed) then the particle gets deflected to the right and Aip 
decreases to a value n > Acpu) > 0. If we now increase p 
to a value p(b) > p c (keeping I fixed) the deflection angle 
Aip decreases to a value < Aip^) < A</?( a ) . If we then 
decrease I (keeping p fixed) the deflection angle decreases 
further. At I — the particle comes in along the y-axis, 
slows down, but now has sufficient energy to cross the 
barrier and move towards y = +oo. This orbit has Aip = 
0. If I is decreased further to a negative value Z( c ) < 
the particle gets deflected to the left giving a negative 
deflection angle —it < A</?( c ) < 0. If p is then decreased 
to a value pt^ < p c (keeping I at Im < 0) the defection 
angle decreases further to a value — it < AipM) < A</?( c ). 
If I is increased (keeping p fixed at P(d)) the deflection 
angle decreases even further and at I = it reaches the 
value Aip = —tt. Upon returning to the starting point of 
closed path 7 in the (Z,p)-plane that encircles the critical 
point (l,p) = (0,p c ) the deflection angle thus is increased 
by27r, 



dAW(l,p) dAW(l,p) 



01 



-dl 



dp 



dp = 2n. (14) 



Such a non-zero value from a closed loop 7 only occurs 
when the critical point (0,p c ) is encircled by 7. If a 
system has loops of regular values for which the deflection 
angle is increased by (multiples of) 2ir we say the system 
has scattering monodromy. 

For the classical system, the angular momentum L and 
the Hamiltonian function H are two independent con- 
stants of motion which are in involution. The classical 



system is therefore integrable. The level set {L = l,H = 
E} in phase space for a regular value (I, E) topologically 
is a cylinder. The invariant cylinder K x S 1 consists of 
an orbit (~ M) as shown in Fig. [5] and all its partners 
with different angle of incidence (~ S* 1 ) but the same 
I and E. At the critical value (0, E c ) the gradients of 
L and H are linearly dependent. The critcial level set 
{L = 0.H = E c } is topologically a cone. It consists 
of the equilibrium point at the origin and all orbits ap- 
proaching it forward or backward in time. The invariant 
cone (a pinched cylinder) in the phase space of a scat- 
tering system is the analogue of the pinched torus in a 
bound system. A loop of invariant cylinders that encir- 
cles the invariant cone cannot be contracted, and as a 
result Aip shows scattering monodromy. 

Conclusions. — We have shown that the quantum scat- 
tering phase shift S for a smooth radially symmetric re- 
pulsive potential cannot be globally defined. The clas- 
sical analogue is that the deflection angle changes by 
2-7T upon traversing a loop in the space of constants of 
motion that encloses the critical value corresponding to 
the equilibrium point. As opposed to the more ab- 
stract consequences of monodromy in compact systems 
this phenomenon is 'directly' observable, e.g. by play- 
ing marbles (neglecting moments of inertia) on a surface 
with a rotationally symmetric bump. 



[4. 
[5 

[6 
[7; 

[8. 
[9 

[10 

[11 



H. Goldstein, Classical Mechanics (Addison-Wesley, 
Reading, MA, 1980), 2nd ed. 
A. Einstein, Verh. DPG 19, 82 (1917). 
V. I. Arnold, Mathematical Methods of Classical Mechan- 
ics, vol. 60 of Graduate Texts in Mathematics (Springer, 
Berlin, 1978). 

J. J. Duistermaat, Comm. Pure Appl. Math. 33, 687 
(1980). 

R. H. Cushman and L. M. Bates, Global Aspects of 
Classical Integrable Systems (Birkhauser, Basel, Boston, 
Berlin, 1997). 

R. H. Cushman and J. J. Duistermaat, Bull. Amer. Math. 
Soc. 19, 475 (1988). 

R. H. Cushman and D. A. Sadovskii, Europhys. Lett. 47, 
1 (1999). 

H. Waalkens, H. R. Dullin, and P. H. Richter, Physica D 
196, 265 (2004). 

R. H. Cushman, H. R. Dullin, A. Giacobbe, D. D. Holm, 
M. Joyeux, P. Lynch, D. A. Sadovskii, and B. Zhilinskii, 
Phys. Rev. Lett. 93, 024302 (2004). 

D. A. Sadovskii and B. I. Zhilinskii, Phys. Lett. A 256, 
235 (1999). 

N. Froman and P. Froman, JWKB Approximation: 
Contributions to the Theory (North-Holland Publishing 
Company, Amsterdam, 1965). 



